Applied  Thermal  Engineering  76  (2015)  64-75 


ELSEVIER 


Contents  lists  available  at  ScienceDirect 

Applied  Thermal  Engineering 

journal  homepage:  www.elsevier.com/locate/apthermeng 


Applied 

Thermal 

Engineering 


Design  .  Processes  .  Equipment  -  Economics 


Research  paper 

Effects  of  furnace  chamber  shape  on  the  MILD  combustion 
of  natural  gas 

Yaojie  Tu,  Hao  Liu  ,  Sheng  Chen,  Zhaohui  Liu,  Haibo  Zhao,  Chuguang  Zheng 

State  Key  Laboratory  of  Coal  Combustion,  Huazhong  University  of  Science  and  Technology,  Wuhan  430074,  PR  China 


CrossMark 


HIGHLIGHTS 


•  The  impact  of  furnace  body  configuration  on  MILD  combustion  was  studied. 

•  Stronger  recirculation  and  smaller  temperature  gradient  are  achieved  at  larger  divergence  angle. 

•  Larger  backflow  and  chemical  flame  occupations  are  obtained  at  larger  divergence  angle. 

•  NO  emission  reduces  at  larger  divergence  angle. 

•  Larger  divergence  angle  is  helpful  for  establishing  MILD  combustion. 
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The  objective  of  this  work  is  to  assess  the  effects  of  furnace  chamber  shape  (by  changing  the  angle  a 
between  the  furnace  roof  and  sidewall)  on  the  MILD  (moderate  or  intense  low-oxygen  dilution)  com¬ 
bustion  characteristics  for  natural  gas  using  high  temperature  air  (1573  K)  with  the  aid  of  CFD  simula¬ 
tions.  The  present  simulations  agree  well  with  the  previous  experimental  results,  and  some  new  findings 
are  revealed.  For  instance,  when  increasing  angle  a  under  the  present  furnace-burner  conditions, 
stronger  recirculation,  smaller  temperature  gradient  as  well  as  larger  reaction  zone  are  obtained,  indi¬ 
cating  the  easier  establishment  of  MILD  combustion.  In  contrast,  the  combustion  regime  transforms  from 
MILD  combustion  towards  flame  combustion  in  the  stability  limits  of  diffusion  combustion  provided  by 
Wiinning  et  al.  (Prog.  Energy  Combust.  Sci.  23  (1997)  81-94)  when  reducing  angle  a.  It  is  also  suggested 
that  there  is  a  critical  value  for  angle  a  below  which  MILD  combustion  cannot  be  established  under  the 
present  conditions. 

©  2014  Elsevier  Ltd.  All  rights  reserved. 


1.  Introduction 

Moderate  or  intense  low-oxygen  dilution  (MILD)  combustion  [1  ] 
has  become  one  of  the  most  promising  combustion  technologies  in 
industrial  heating  for  low  NOx  emission  and  high  efficiency  of 
combustion.  MILD  combustion  is  featured  as  intense  recirculation 
of  hot  flue  gas  inside  the  furnace,  which  on  one  hand  improves  the 
temperature  of  reactants  above  their  ignition  point,  on  the  other 
hand  dilutes  the  mixture  of  reactants  before  combustion.  In  this 
regard,  combustion  takes  place  in  an  environment  of  highly  diluted 
reactants  where  auto-ignition  is  also  allowed  to  happen.  The 
combustion  process  under  MILD  condition  is  very  different  from 
conventional  combustion  since  the  internal  recirculation  is  so 
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strong  that  the  combustion  reaction  zone  spreads  to  a  broader 
space.  Hence,  both  the  temperature  and  the  chemical  species  are 
distributed  homogeneously,  and  the  peak  temperature  as  well  as 
NOx  emission  is  reduced  as  well. 

This  technology  is  also  named  as  high  temperature  air  com¬ 
bustion  (HiTAC)  when  highly  preheated  temperature  air  is  used 
(usually  >1000  °C)  [2,3],  or  flameless  combustion  or  flameless 
oxidation  (FLOX)  [4,5],  or  colorless  distributed  combustion  (CDC) 
[6,7  .  For  convenient  description,  the  term  “MILD  combustion”  is 
used  uniformly  throughout  the  present  paper.  Plenty  of  valuable 
work  has  been  done  on  MILD  combustion  by  such  as  Cavaliere  et  al. 
[1  ,  Tsuji  et  al.  [2,3],  Wunning  et  al.  [4,5  ,  Gupta  et  al.  6,7  ,  Webber 
et  al.  [8-13],  Blasiak  et  al.  14-16],  Dally  et  al.  17,18]  and  Mi  et  al. 
[19-23]  Moreover,  detailed  descriptions  for  MILD  combustion  and 
HiTAC  have  been  respectively  reviewed  by  Cavaliere  et  al.  1  ]  and 
Tsuji  et  al.  [3  To  summarize,  Wunning  et  al.  gave  the  major  clari¬ 
fications  of  MILD  combustion  in  ref  5],  and  demonstrated  that  the 
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Nomenclature 

A- 

reference  furnace  cross  sectional  area 

Aa 

absolute  furnace  cross  sectional  area 

a 

angle  between  furnace  roof  and  sidewall 

me 

entrained  flue  gas  mass  flow  rate 

Teh 

characteristic  time  scale  of  chemical  reaction 

mf 

fuel  mass  flow  rate 

^turb 

characteristic  time  scale  of  turbulent  mixing 

ma 

air  mass  flow  rate 

Yf 

mass  fraction  of  species  (i) 

* 

equivalent  internal  recirculation  rate 

I<v 

V 

velocity  vector 

Ro 

oxidation  mixture  ratio 

ji 

mass  diffusion  flux  of  species  (i) 

m0 

local  oxygen  mass  fraction 

Ri 

net  production  rate  of  species  (i) 

m'o 

oxygen  mass  fraction  needed  to  complete  combustion 

Vx 

x  component  of  mean  velocity 

Rc 

critical  oxidation  mixture  ratio 

I<v 

P 

internal  recirculation  rate 
mixture  density 

Rf 

chemical  flame  occupation  degree 

realization  of  MILD  combustion  has  nothing  to  do  with  operational 
conditions  like  fuel  types,  air  or  fuel  preheating  temperature,  nor 
burner  or  furnace  configurations,  as  long  as  there  exits  an  intense 
recirculation  of  hot  combustion  products  inside  the  furnace.  So  far, 
MILD  combustion  has  been  applied  in  the  industrial  furnaces, 
especially  in  steel  making.  However,  R&D  of  MILD  combustion 
mainly  deals  with  gaseous  fuels  and  light  oils  due  to  their  prefer¬ 
able  combustion  characteristics  to  solid  fuels,  but  the  latest  studies 
of  MILD  combustion  with  solid  fuels  like  pulverized  coal  and 
biomass  still  can  be  found  elsewhere  [8,11-13,16,18,24-36  . 

It  is  worth  mentioning  that  the  IFRF  (International  Flame 
Research  Foundation)  [8-12]  has  carried  out  a  series  of  MILD 
combustion  experiments  for  natural  gas,  light  oil,  heavy  oil  and 
pulverized  coal  with  high  temperature  air,  which  leads  to  a  deep 


insight  into  the  fundamental  and  practical  issues  of  MILD  com¬ 
bustion.  All  the  tests  by  IFRF  were  conducted  on  a  0.58  MW  vertical 
furnace  with  the  burner  inlets  and  the  exhaust  gas  exit  located  at 
the  opposite  side  of  the  furnace  walls.  Later,  a  series  of  MILD 
combustion  experiments  with  natural  gas,  sawdust  and  pulverized 
coal  under  both  O2/N2  and  O2/CO2  atmospheres  were  carried  out  at 
Adelaide  [17,18]  on  a  20  KW  recuperative  furnace  with  the  burner 
inlets  and  the  exhaust  gas  exit  located  at  the  same  side  of  the 
furnace  walls.  According  to  the  experimental  and  simulation  results 
of  the  above  two  furnaces  which  have  different  geometrical  con¬ 
figurations,  the  internal  recirculation  is  usually  found  to  be  stronger 
inside  the  Adelaide's  furnace  as  compared  to  the  IFRF's  furnace.  The 
reason  behind  this  is  that  for  Adelaide's  furnace,  the  combustion 
products  would  totally  flow  back  into  the  upstream  furnace;  while 


100 


(b)  (c) 


Fig.  1.  Schematic  diagram  of  the  reference  furnace-burner  system:  (a)  the  furnace;  (b)  the  burner;  (c)  top  view  of  the  furnace  (unit:  mm). 


66 


Y.  Tu  et  al.  /  Applied  Thermal  Engineering  76  (2015)  64—75 


Table  1 

Experimental  conditions  in  IFRF  test  9  . 


Flow  rate 
(kg/h) 

Temperature 

(I<) 

Enthalpy 

(MW) 

Components 

(vol.%) 

Fuel 

47 

298 

0.58 

CH4  =  87.8,  C2H6  =  4.6, 
C3H8  =  1.6,  N2  =  5.5 

Oxidizer 

830 

1573 

0.35 

02  =  19.5,  N2  =  59.1, 

H20  =  15,  C02  =  6.4, 

NO  =  110  ppm 

Exhaust 

gases 

877 

1493 

0.38 

02  =  1.6,  NO  =  140  ppm 

for  IFRF's  furnace,  only  part  of  the  combustion  products  recedes 
backward.  This  to  some  extent  indicates  that  the  internal  configu¬ 
ration  of  furnace  or  its  exit  location  can  be  an  important  factor  on 
the  internal  dynamic  flow  field,  but  not  only  the  burner  system. 

Over  the  recent  years,  the  investigations  on  MILD  combustion 
which  were  related  to  the  geometrical  parameters  mainly  focused 
on  the  burner  system  17-22,35-38]  (e.g.  injection  velocity  or  di¬ 
rection,  nozzle  separation  distance,  mixing  pattern,  burner  location 
and  arrangement).  Both  the  experiments  and  numerical  simula¬ 
tions  were  carried  out  to  study  the  above  parameters.  However,  less 
work  has  been  done  with  the  furnace  internal  configuration, 
namely  the  combustion  chamber  shape.  On  this  side,  Schaffel  N. 
et  al.  13]  proposed  a  conceptual  design  of  a  supercritical  boiler 
utilizing  MILD  combustion  technology,  and  they  found  that  the 
varied  boiler  body  shapes  made  a  difference  to  the  internal  recir¬ 
culation,  coal  burnout  ratio  as  well  as  the  temperature  distribution 
(which  is  a  dominant  factor  on  thermal  NOx  formation).  By  properly 
designing  the  boiler  parameters,  higher  combustion  efficiency  and 
heat  flux  as  well  as  lower  NOx  emission  could  be  realizable  under 
MILD  combustion  in  comparison  with  conventional  boiler.  In  this 
regard,  the  geometrical  parameters  of  the  burner  as  well  as  the 
combustion  chamber  itself  can  play  very  important  role  on  the 
furnace  internal  recirculation,  which  is  the  essential  condition  for 
achieving  MILD  combustion. 

From  the  academic  point  of  view,  flame  structure  can  be  influ¬ 
enced  or  even  dominated  by  the  dynamic  flow  inside  the  furnace. 
However,  the  internal  dynamic  flow  is  significantly  depended  on 
the  characteristics  of  both  burner  and  combustion  chamber  con¬ 
figurations.  Up  to  now,  there  are  kinds  of  burners  designed  for  or¬ 
dinary  or  specific  applications,  i.e.  direct  flow  burner,  swirl  burner, 
staged  burner,  and  MILD  burner.  The  different  design  concepts  of 
the  burners  rely  on  the  different  demands,  for  example,  swirl 
burner  is  featured  as  central  recirculation  of  hot  flue  gas  to  enhance 
the  ignition  stability,  while  MILD  burner  is  applied  to  obtain  a 
uniform  temperature  field  and  lower  NOx  emission.  As  for  furnace 


Table  2 

Furnace  geometrical  parameters  of  the  considered  cases. 


Casel 

Case2 

Case3 

Case4 

Case5 

a 

o 

o 

00 

85° 

o 

o 

CD 

95° 

100° 

Li  (mm) 

2998 

2524 

2000 

1426 

795 

L2  (mm) 

795 

1426 

2000 

2524 

2998 

chamber  configuration,  normal  patterns  of  cylinder  or  cubic  shapes 
are  mostly  used  for  researches.  The  regime  from  conventional 
combustion  to  MILD  combustion  is  transformed  by  simply  replac¬ 
ing  the  corresponding  burners.  However,  the  potential  impact  of 
furnace  chamber  shape  on  combustion  characteristics  is  barely 
studied,  especially  under  MILD  operation.  On  the  other  hand, 
reheating  furnace  is  one  of  the  top  energy  consuming  industrial 
devices,  thus  it  is  necessary  to  optimize  the  geometric  parameters 
to  obtain  better  performance  and  save  energy,  especially  when  the 
function  of  a  furnace  changes  or  new  furnace  is  going  to  be  built. 

The  main  objective  of  the  present  study  is  to  investigate  the 
effects  of  furnace  chamber  shape  (by  changing  the  angle  a  between 
the  furnace  roof  and  sidewall)  on  MILD  combustion  with  the  aid  of 
CFD  simulation  based  on  the  IFRF  experimental  instrument.  First, 
experimental  data  from  the  reference  case  (a  =  90°)  was  employed 
to  validate  the  present  simulations.  Then,  the  characteristics  of 
internal  mixing,  combustion  as  well  as  NOx  emission  were 
respectively  compared  for  the  considered  cases.  This  study  has 
figured  out  the  change  rules  of  the  in-furnace  flow  fields  and 
combustion  behaviors  against  the  angle  a ,  and  can  provide  some 
useful  information  on  optimizing  furnace  geometric  parameters 
under  MILD  operation. 

2.  Computational  details 

2.2.  Reference  conditions  and  considered  furnace  shape 

The  MILD  combustion  experiment  with  natural  gas  of  the 
reference  case  were  performed  by  IFRF  in  a  semi-industrial  scale 
furnace  (see  Fig.  1)  with  a  2  m  x  2  m  cross  section  and  a  length  of 
6.25  m  [9,10  .  The  burner  system  was  consisted  of  a  central  oxidizer 
nozzle  and  two  off-axis  fuel  nozzles.  The  internal  diameter  of  the 
oxidizer  and  fuel  nozzles  were  124  mm  and  10  mm,  respectively. 
The  distance  between  the  oxidizer  and  fuel  nozzles  was  280  mm. 
The  furnace  was  operated  at  0.58  MW  fuel  input  with  the  oxidizer 
preheated  to  1573  K.  Table  1  lists  the  experimental  conditions  of  the 
IFRF  test.  Note  that,  the  high  temperature  oxidizer  was  generated 
from  a  pre-combustor,  therefore,  a  small  quantity  of  H20  and  C02 


Fig.  2.  Schematic  diagram  of  the  considered  furnace  chamber  shape. 


Y.  Tu  et  al  /  Applied  Thermal  Engineering  76  (2015)  64—75 


67 


Table  3 

Global  reaction  mechanism  of  NG  used  in  the  present  study  (units  in  SI). 


Index 

Reactions 

A  [kmol/(m3  s  1)] 

Ea  (J/kmol) 

Reaction  orders 

R.l 

CH4  +  0.5O2  - 

->  CO  +  2H2 

4.4  x  1011 

1.255  x  108 

[CH4]°'7[02]0'8 

R.2 

CO  +  0.502  -> 

co2 

2.24  x  1012 

1.703  x  108 

[CO][O2]025[H2O]05 

R.3 

H2  +  0.502  -> 

h2o 

5.69  xlO11 

1.46  x  108 

[H2][02][H20] 

R.4 

C2H6  +  3.502 

2C02  +  3H20 

6.186  x  109 

1.256  x  108 

[c2h6]01[o2]165 

R.5 

C3H8  +  502  —5 

►  3C02  +  4H20 

6.186  x  109 

1.256  x  108 

[C3H8]01[O2]1’65 

a  The  reaction  rate  coefficient  k  =  Aexp(-EaIRuT),  where  Ru  =  8315  J  kmol  1  I( 


Fig.  3.  Computational  grid  for  simulating  the  reference  case  with  565224  cells. 


were  contained  in  the  oxidizer  together  with  a  relatively  large 
quantity  of  NO. 

The  considered  furnace  chamber  shapes  in  the  current  study 
were  varied  by  changing  the  angle  a  between  the  furnace  roof  and 
sidewall  as  illustrated  in  Fig.  2.  Five  angles  were  examined 
including  two  acute  angles  (a  =  80°,  and  a  =  85°),  a  right  angle 
(a  =  90°,  the  reference  case),  and  two  obtuse  angles  (a  =  95°,  and 
a  =  100°).  Note  that,  the  total  length  and  the  volume  of  the  furnace 
were  kept  unchanged  in  order  to  have  a  same  average  volume  heat 


load  as  well  as  residence  time.  The  side  length  of  the  furnace  roof 
and  floor  were  summarized  in  Table  2  together  with  the  corre¬ 
sponding  angle. 

2.2.  Simulation  approaches 

The  simulations  were  performed  with  the  aid  of  the  commercial 
CFD  (computational  fluid  dynamics)  software  ANSYS  Fluent  14.0. 
The  viscous  turbulent  flow  was  solved  using  standard  k-s  model 
and  standard  wall  function  [39  .  The  PI  model  with  the  weighted 
sum  of  gray  gas  model  (WSGGM-domain-based)  was  used  to  solve 
the  radiation  and  absorption  between  the  gaseous  molecules  and 
furnace  walls  [24,40].  The  SIMPLE  scheme  was  used  for  pressur¬ 
e-velocity  coupling  [41  .  The  main  challenge  in  simulating  MILD 
combustion  is  to  select  a  proper  combustion  model.  Since  the  mean 
reaction  rate  under  MILD  condition  is  slower  than  conventional 
combustion,  the  characteristic  time  scale  of  chemical  reaction  (rCh) 
becomes  comparable  to  that  of  turbulent  mixing  (rturb),  thus  the 
turbulence-chemistry  interactions  should  be  treated  with  finite 
rate  approaches  [42].  In  this  regard,  the  eddy-dissipation  model 
(EDM)  and  the  eddy-dissipation  concept  (EDC)  are  expected  to  be 
considered  in  modeling  MILD  combustion.  However,  in  the  present 
CFD  codes,  the  EDM  allows  only  a  two-step  global  mechanism, 
while  the  EDC  allows  both  global  and  detailed  mechanisms. 
Moreover,  M.  Vascellari  et  al.  [24  have  demonstrated  that  the 
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0 
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Fig.  4.  Grid-independence  test  in  axial  temperature,  velocity,  02  and  C02  distributions  for  the  reference  case. 
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coupling  of  EDC  with  detailed  mechanism  would  produce  better 
results  of  the  MILD  combustion.  Therefore,  we  choose  EDC  as  the 
combustion  model  in  the  present  study  rather  than  EDM. 

In  the  EDC  model,  the  transporting  equations  for  the  chemical 
species  are  governed  by: 

^(pYO  +  VjWYj)  =-V-X  +Rt  (1) 

where  p  is  the  mixture  density,  Yj  is  thejocal  mass  fraction  of  each 
species  (i),  ~v  is  the  velocity  vector,  and  Ji  is  the  mass  diffusion  flux 
of  species  (i),  and  Ri  is  the  net  rate  of  production  of  species  (i)  by 
chemical  reaction.  The  mass  diffusion  flux  Ji  and  the  chemistry 
source  term  Ri  of  ith  species  are  calculated  by: 


where  Sct  is  the  turbulence  Schmidt  number,  DI>m  is  the  mass 
diffusion  coefficient  for  each  species,  £  is  the  size  of  fine  scales,  and 
Y*  is  the  fine-scale  species  mass  fraction  after  reacting. 

In  order  to  reduce  the  calculation  consumption,  a  five-step 
global  reaction  mechanism  is  used  in  the  present  study  for  the 
fuel  as  listed  in  Table  3.  The  global  mechanism  contains  three  re¬ 
actions  for  CH4  and  one  reaction  for  C2H6  and  C3H8,  respectively.  To 
examine  the  suitability  of  the  global  reaction  mechanism  on 
modeling  MILD  combustion,  the  simulations  of  the  reference  case 
with  a  detailed  chemical  kinetic  mechanism  for  CH4  from  GRI- 
Mech  3.0  [43]  is  also  carried  out.  As  indicated  by  M.  Vascellari 
et  al.  [24  ,  the  low  NOx  concentration  has  a  reduced  impact  on  the 
other  flow  properties,  like  temperature  and  species  distributions  as 
well  as  mass,  momentum  and  energy  balances.  Moreover,  due  to 
the  longer  characteristic  time  scale  of  NOx  kinetic  with  respect  to 
the  combustion  kinetics,  the  NOx  production  can  be  evaluated  by 
post-processing  approach  24].  Finally,  the  N-containing  reactions 
are  removed  in  GRI-Mech  3.0  mechanism,  and  only  219  reversible 
reactions  with  36  chemical  species  are  calculated. 

The  fuel,  oxidizer  and  flue  gas  boundary  conditions  were  set  as 
constant  mass  flow  inlet  and  pressure  outlet,  respectively.  The 
furnace  wall  was  set  at  constant  temperature  of  1500  K  with  a 
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Fig.  5.  Comparison  of  the  experimental  A )  and  simulated( - global  mechanism, - GRI  3.0  mechanism)  results  of  the  reference  case:  (a)  velocity:  (b)  temperature;  (c)  oxygen 

concentration:  (d)  methane  concentration  at  x  =  0.15  m,  0.73  m,  2.05  m,  z  =  0  m. 
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constant  internal  emissivity  rate  of  0.6.  Convergence  was  obtained 
when  the  residence  of  continuity  was  lower  than  10”5,  and  10~6  for 
the  other  variables  (velocity,  energy,  PI  and  species,  etc.).  In  addi¬ 
tion,  the  maximum  temperature  and  velocity  on  the  symmetrical 
plane  (z  =  0)  was  monitored,  regarded  as  another  criterion  for 
convergence  when  their  variations  were  less  than  1  I<  and  0.1  m/s, 
respectively. 

3.  Results  and  discussion 

3  A.  Grid-independent  test 

Owning  to  the  symmetrical  geometry  of  the  present  furnace, 
only  one  quarter  of  the  furnace  was  modeled.  Fig.  3  shows  the 


computational  grid  used  for  simulating  the  reference  case,  which  is 
consisted  of  565224  hexahedral  cells  with  size  functions  on  the 
furnace  roof  and  sidewall.  In  order  to  check  the  independence  of 
the  grid,  simulation  work  was  also  conducted  for  the  reference  case 
with  a  fine  grid  consisted  of  1319902  cells.  Refinement  on  the 
furnace  roof  and  sidewall  was  made  to  the  grid.  Fig.  4  compares  the 
gas  temperature,  velocity,  O2  and  CO2  distribution  profiles  along  the 
furnace  centerline  from  the  two  grids.  As  demonstrated,  no  obvious 
difference  in  the  axial  temperature,  velocity  and  species  distribu¬ 
tions  can  be  distinguished  between  the  coarse  grid  and  the  fine 
grid.  In  order  to  save  the  computational  time,  the  coarse  grid  with 
565224  cells  was  utilized  for  simulating  the  reference  case,  and 
similar  meshing  method  and  grid  size  were  employed  for  the 
subsequent  cases. 


Fig.  6.  Contours  of  x  component  of  mean  velocity  (a):  x-y  plane  (z  =  0);  (b):  x  =  1000  mm  plane;  (c):  x  =  2000  mm  plane.  The  black  curves  denote  the  level  of  Vx  =  0. 
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Fig.  7.  Schematic  diagram  of  the  velocity  field  (up)  and  aerodynamic  flow  field  (down) 
inside  the  furnace  for  the  reference  case. 


3.2.  Validation  of  the  computational  models 

In  order  to  examine  the  adequacy  of  the  adopted  models  with 
the  global  reaction  mechanism  in  the  present  study,  numerical 
simulations  were  also  carried  out  using  the  detailed  mechanism 
from  GRI-Mech  3.0  for  the  reference  case.  The  simulated  results  of 
mean  velocity,  temperature,  02  and  CH4  concentrations  in  furnace 
from  both  global  and  detailed  mechanisms  were  compared  with 
the  experimental  data  [9,10]  at  different  measurement  locations 
(see  Panels  a-d  of  Fig.  5). 

As  displayed  in  Fig.  5,  the  calculated  results  from  both  global  and 
detailed  mechanisms  have  shown  good  agreement  with  the 
measured  data  for  all  of  the  above  parameters.  However,  the 
simulated  results  of  temperature  and  CH4  concentration  distribu¬ 
tions  from  the  global  mechanism  seem  to  be  unsatisfactory  near 
the  fuel  injection  zone  (Y  =  280  mm).  In  addition,  the  02  con¬ 
sumption  is  found  to  be  improved  with  the  global  mechanism.  The 
reason  behind  these  differences  is  probably  that,  the  simplified 
global  mechanism  has  not  taken  the  intermediate  species  (i.e.  CH, 
OH,  etc.)  into  account,  and  local  details  of  the  chemical  reaction 
cannot  be  precisely  captured  especially  in  the  fuel  rich  region. 
Nevertheless,  the  characteristics  of  the  dynamic  flow  field  and 
combustion  status  still  can  be  generally  expressed  with  the  global 
mechanism.  Since  the  computing  cost  increases  remarkably  with 
the  GRI-Mech  3.0  mechanism,  the  simplified  global  mechanism  is 
adopted  to  save  the  computing  consumption. 


Fig.  8.  Furnace  temperature  against  I<v  for  the  present  simulation  cases  and  the 
schematic  diagram  of  stability  limits  for  methane  in  diffusion  combustion  from  Ref.  [4]. 


3.3.  Effects  of  furnace  chamber  shape  on  mixing 

MILD  combustion  is  characterized  as  intense  recirculation  of  flue 
gas  inside  the  furnace,  and  Fig.  6  shows  the  contours  of  the  x 
component  of  mean  velocity  (\/x)  in  the  x-y  plane  (z  =  0)  and  the  y-z 
planes  at  x  =  1000  mm  and  2000  mm  for  the  five  cases.  The  border 
(\4  =  0)  between  the  upward  and  downward  velocity  can  be 
distinguished  by  the  black  curves  and  the  contour  legend.  Similar 
flow  field  distributions  can  be  observed  for  the  five  cases.  Concretely, 
they  all  have  a  central  upward  flow  ( Vx  >  0)  and  downward  back- 
flows  (Vx  <  0)  near  the  furnace  sidewalls  and  four  corners.  Note  that, 
the  backflow  even  can  be  observed  in  the  furnace  central  region 
when  a  =  80°.  As  observed,  with  an  increase  of  angle  a ,  the  position 
where  backflow  starts  occurring  in  x-y  plane  (z  =  0)  is  moving  to¬ 
wards  upstream  furnace.  Meanwhile,  the  occupation  of  backflow  in 
x-y  plane  is  becoming  larger  as  angle  a  increases.  Especially  when 
a  =  100°,  the  backflow  occupies  the  whole  surrounding  area  of  the 
central  upward  flow  throughout  almost  the  whole  furnace. 

Fig.  7  illustrates  the  simplified  schematic  diagram  of  the  in¬ 
furnace  flow  and  aerodynamic  fields.  In  order  to  identify  the  dilu¬ 
tion  degree  of  reactants  by  entrained  flue  gas,  here,  the  internal 
recirculation  rate  I<y  is  introduced,  which  is  first  defined  by  Wiin- 
ning  et  al.  [4]  as  the  ratio  of  mass  flow  rate  of  the  entrained  flue  gas 
(me)  to  that  of  the  initial  fuel  (m/)  and  air  (ma)  jets  at  a  certain  plane 
(see  Fig.  7).  It  is  expressed  as: 
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Fig.  9.  Contours  of  temperature  fields  in  the  x-y  plane  (z  =  0)  for  the  five  cases. 
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ff  pvxdydz  -ma-mf 
K  Trie  mup  —  rna  —  ttij-  A(up) 

v  ~  (ma  +  xrif)  ~  (ma  +  mf)  ~  (ma  +  mf) 

(4) 

where  A(up)  is  the  total  area  in  a  certain  cross  section  with  an 
upward  axial  velocity  (vx>0),  and  is  represented  by  the  central  part 
within  the  black  curves  as  shown  in  Fig.  6b— c.  Note  that,  the  cross 
sectional  area  at  a  fixed  axial  distance  varies  with  different  angle  a , 
particularly,  in  the  upstream  furnace  (x  <  3125  mm),  the  cross 
sectional  area  reduces  as  angle  a  increases.  To  eliminate  this 

geometrical  discrepancy,  the  equivalent  internal  recirculation  rate 

* 

Kv  is  used  here,  and  is  expressed  as: 

*  Ar 

Kv  =  ~r~  Kv  (5) 


x  =  1000  mm  and  x  =  2000  mm  planes  versus  the  furnace  tem¬ 
perature  for  the  five  cases  is  also  plotted  in  Fig.  8. 

From  Fig.  8,  it  can  be  found  that  the  main  reaction  zone  (see 
Fig.  12)  is  generally  in  the  regime  of  MILD  combustion  for  all  of  the 
five  cases,  except  for  easel  and  2  in  x  =  1000  mm  plane.  Hence,  it  is 
revealed  that  there  exists  a  critical  value  for  angle  a  below  which 
the  steady  MILD  combustion  would  not  occur  under  the  present 
conditions.  Furthermore,  it  implies  that  when  establishing  MILD 
combustion  in  a  furnace  with  a  relatively  small  angle  a ,  the  initial 
injection  conditions  might  need  to  be  intensified  (such  as 
increasing  injecting  velocity  or  preheating),  while  in  a  furnace  with 
a  relatively  large  a ,  the  initial  injection  condition  could  be  relaxed, 
and  it  is  also  helpful  for  achieving  MILD  combustion. 

In  addition,  the  combustion  regime  movement  against  angle  a  is 
observed  of  a  similar  rule  in  both  x  =  1000  mm  and  x  =  2000  mm 

planes  according  to  Fig.  8.  Specifically,  the  larger  a  always  leads  to 

* 

lower  temperature  and  improved  Kv  at  a  certain  cross  section,  this 
trend  can  be  clearer  when  a  =  100°.  It  indicates  that  increasing 
angle  a  can  strengthen  the  mixing  between  the  reactants  and  flue 
gas.  This  is  resulted  from  the  reduced  cross  sectional  area  with  the 

larger  angle  a  which  promotes  the  mean  velocity  in  the  upstream 

* 

furnace.  Also  larger  I<v  is  observed  in  x  =  2000  plane  than 
x  =  1000  mm  plane  with  the  same  furnace  shape.  According  to  Eq. 
(4),  it  can  be  seen  that  I<v  is  in  direct  proportion  to  the  product  of 
mean  density,  velocity  and  area  of  the  upward  flow.  However,  the 
variance  of  mean  density  is  negligible  compared  to  that  of  mean 
velocity,  and  the  occupation  ratio  of  the  upward  flow  is  found  to 
decrease  as  angle  a  increases  according  to  Fig.  6b-c  (normalized  by 
the  reference  case).  Hence,  it  implies  that  the  increase  of  mean 
velocity  is  the  dominant  factor  on  the  strengthened  internal  recir¬ 
culation.  It  further  confirms  that  when  designing  a  furnace  under 
MILD  operation,  a  larger  angle  a  would  be  helpful  for  establishing 
MILD  combustion. 


where  Aa  is  the  absolute  cross  sectional  area  at  a  certain  axial  dis¬ 
tance,  and  Ar  is  the  reference  cross  sectional  area  at  the  same  axial 
distance.  In  the  present  study,  case3  is  the  reference  case  with  a 
constant  cross  sectional  shape  along  the  furnace  centerline.  Thus, 

Ar  =  lm2  (one  quarter  of  the  geometry  is  simulated)  is  used  for  the 

* 

Kv  calculations. 

Fig.  8  shows  the  relationship  between  the  I<y  and  the  furnace 
temperature  in  diffusion  combustion  for  methane,  which  was 
provided  by  Wiinning  et  al.  4].  According  to  Fig.  8,  the  flame  sta¬ 
bility  limit  is  consisted  of  three  regimes  (regions  A,  B,  C  stands  for 
stable  conventional  combustion,  unsteady  transitional  regime, 
MILD  combustion).  Since  methane  takes  up  the  largest  proportion 
(87.8%  by  volume)  in  the  present  fuel,  the  relationship  in  Fig.  8 
could  be  applied  in  order  to  identify  the  combustion  regime  of 

the  present  cases  regardless  of  slight  difference  in  fuel  composition. 

* 

Therefore,  the  corresponding  equivalent  recirculation  rate  Kv  in 


Fig.  11.  Schematic  diagram  of  chemical  flame. 


3.4.  Effects  of  furnace  chamber  shape  on  combustion 

The  uniform  temperature  distribution  has  impelled  the  appli¬ 
cation  of  MILD  combustion  in  the  industrial  heating.  The  intuitive 
way  to  inspect  the  uniformity  of  the  temperature  distribution  is  to 
obtain  the  temperature  field  inside  the  furnace  directly.  According 
to  the  experiments  [11],  the  combustion  temperature  peak  did  not 
exceed  1550  °C  (1823  K),  and  the  whole  furnace  was  filled  with 
combustion  products  with  temperature  in  the  range  of 
1300-1500  °C  (1573-1773  K).  Fig.  9  displays  the  temperature  fields 
in  the  x-y  plane  (z  =  0)  for  the  five  cases.  The  isotherms  of  1600  K, 
1700  I<  and  1800  I<  are  depicted  in  the  temperature  fields.  A  black 
arrowed  line  is  also  given  to  represent  the  variation  trend  of  the 
area  with  temperature  above  1600  K.  It  can  be  found  that,  with  the 
expansion  of  angle  a ,  the  axial  distance  of  the  isotherm  of  1600  K 
reaches  further.  Moreover,  the  maximum  gas  temperature  inside 
the  furnace  is  shown  in  Fig.  10  for  the  five  cases,  and  the  temper¬ 
ature  peak  is  observed  to  reduce  with  the  increase  of  angle  a.  In  this 
regard,  it  can  be  known  that  the  temperature  gradient  is  lower  and 
the  temperature  uniformity  improves  as  angle  a  increases. 

In  the  IFRF's  experiments,  there  were  no  visible  flames  being 
observed  when  natural  gas  was  fired  ill  ,  thus,  the  conventional 
identification  of  the  flame  size  with  naked  eye  is  no  longer  appli¬ 
cative  to  describe  MILD  combustion.  In  some  sense,  MILD  com¬ 
bustion  is  more  like  the  “diffusive  combustion”  or  “volumetric 
combustion”,  and  the  chemical  reaction  of  the  combustion  re¬ 
actants  takes  place  in  a  broader  space  than  conventional  combus¬ 
tion.  In  order  to  identify  the  volume  of  reaction  zone,  a  concept  of 
chemical  flame  [15,16]  is  introduced  in  the  present  work.  It  is 
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described  by  means  of  the  oxidation  mixture  ratio  ( Ro),  which  is 
expressed  as: 


m0 

m0  +  m'0 


facilitated,  and  the  chemical  reaction  slows  down.  All  the  above 
observations  of  combustion  temperature  and  chemical  flame 
indicate  that,  the  larger  angle  a  would  lead  to  a  larger  reaction  zone 
as  well  as  slower  reaction  rate,  if  more  uniform  temperature  is 
required,  angle  a  is  suggested  to  be  larger. 


where  mo  is  the  local  mass  fraction  of  oxygen,  and  m0  is  the  total 
mass  fraction  of  oxygen  needed  to  complete  combustion.  In  the 
present  study,  the  combustible  species  contains  CH4,  C2H6,  C3H8, 
CO,  H2.  To  complete  the  combustion,  the  needed  oxygen  amount 
(m'0)  is  calculated  as: 

m'0  =  (2YCh4  +  3.5YC2h6  +  5YC3h8  +  0.5YCO  +  0.5YH2)MO2  (7) 

where  Y;  is  the  volume  fraction  of  the  combustible  species,  and  M02 
is  the  molecular  mass  of  oxygen.  Thus,  Eq.  (6)  can  be  transformed 
as: 


3.5.  Effects  of  furnace  chamber  shape  on  NOx  emission 

During  the  combustion  of  natural  gas  under  air  condition,  there 
are  two  major  NO  formation  mechanisms,  namely  the  thermal-NO 
and  the  prompt-NO.  Thermal-NO  is  formed  by  the  oxidation  of 
molecular  nitrogen,  and  is  favored  under  high  temperature  condi¬ 
tion  [44].  Prompt-NO  is  formed  by  the  oxidation  of  nitrogen  radi¬ 
cals  especially  under  fuel-rich  condition  [45].  Besides,  NO  can  be 
reduced  by  the  intermediate  species  like  CH,  radicals,  and  subse¬ 
quently  changes  to  HCN,  which  is  known  as  the  NO-reburning 
mechanism  46].  However,  HCN  can  be  further  oxidized  to  NO.  In 


Ro  = 


yo2m  o2 


Y02M02  +  (2yCH.  +  35yGHfi  +  5VC,H«  +  O.SVco  +  0.5  YHJM0 


'2°6 


3°8 


Y 


0: 


y02  +  2yCH4  +  3-5yC2H6  +  5yC3H8  +  °-5yco  +  0.5YH2 


For  fuel  and  oxidizer  inlet,  the  value  of  Ro  is  0  and  1,  respectively. 
Ro  =  Rc  is  assumed  to  indicate  the  flame  boundary,  and  chemical 
flame  is  assumed  to  be  inside  the  space  where  0  <  Ro  <  Rc  as  shown 
in  Fig.  11.  Note  that,  it  does  not  mean  that  there  is  actual  flame 
inside  the  chemical  flame  boundary,  since  the  aim  of  introducing 
this  concept  is  to  describe  the  reaction  zone  semi  quantitatively.  In 
order  to  validate  the  exact  value  for  Rc ,  experimental  and  simula¬ 
tion  work  has  been  carried  out  by  Yang  et  al.  [15]  and  Li  et  al.  16], 
and  the  value  of  0.99  was  found  to  be  adequate  to  describe  the 
flame  size.  Therefore,  Rc  =  0.99  is  also  adopted  in  this  paper  to 
represent  the  chemical  flame  border. 

Fig.  12  displays  the  contours  of  chemical  flame  shape  in  the  x-y 
plane  (z  =  0)  and  the  y-z  planes  at  x  =  1000  mm  and  2000  mm  for 
the  five  cases.  The  chemical  flame  border  is  represented  by  a  black 
curve.  Similar  flame  pattern  is  observed  with  a  double  parallel  peak 
shape  for  each  case.  According  to  Fig.  12a,  a  slight  increase  of  the 
chemical  flame  length  (the  axial  maximum  distance  of  the  chemical 
flame  boundary)  is  observed  as  angle  a  increases.  Simultaneously, 
the  occupation  ratio  of  chemical  flame  in  both  x  =  1000  mm  and 
x  =  2000  mm  planes  (normalized  by  the  reference  case)  is  found  to 
be  enlarged  with  the  increase  of  angle  a  (see  Fig.  12b-c). 

In  order  to  specify  the  characteristic  of  volumetric  combustion 
under  MILD  condition,  chemical  flame  occupation  degree  (Rp)  is 
introduced  here,  which  is  defined  as  the  ratio  between  the  chem¬ 
ical  flame  volume  (Vf)  and  the  furnace  volume  ( V»,  and  is  calculated 
as: 

Vf 

RF  =  vf  (9) 

vF 

According  to  the  calculated  Rp  for  the  five  cases  (see  Fig.  13),  Rp 
increases  with  angle  a,  which  means  that  the  chemical  flame  vol¬ 
ume  expands  as  the  angle  a  increases.  This  is  caused  by  the 
improved  mean  velocity  and  strengthened  internal  recirculation 
under  larger  angle  a.  Consequently,  the  dilution  of  the  reactants  is 


the  present  simulations,  NO  was  calculated  using  the  post¬ 
processing  method  considering  the  thermal,  prompt  and  the 
reburning  mechanisms. 

Table  4  reports  the  predicted  NO  emission  for  the  five  calculated 
cases.  Note  that,  the  high  temperature  oxidizer  contains  as  high  as 
110  ppm  NO  which  takes  up  a  large  part  of  the  total  NO  emission. 
According  to  Table  4,  the  net  NO  formation  reduces  when 
increasing  the  angle  a.  This  may  be  resulted  from  two  reasons:  First, 
the  reduced  temperature  restrains  the  formation  of  thermal-NO 
since  the  thermal-NO  is  quite  sensitive  to  the  combustion  tem¬ 
perature.  Second,  the  stronger  internal  recirculation  enhances  the 
dilution  degree  of  the  reactants.  Hence,  a  more  reductive  atmo¬ 
sphere  can  even  be  created  which  promotes  the  NO-reburning 
mechanism.  Furthermore,  the  promoted  reductive  atmosphere 
seems  to  be  helpful  for  the  prompt-NO  generation.  However,  it  is 
found  that  the  amount  of  prompt-NO  is  comparable  to  that  of 
reburning-NO.  Consequently,  the  contribution  to  the  NO  reduction 
mainly  comes  from  the  reduction  of  thermal-NO. 

4.  Conclusion 

This  study  numerically  investigated  the  effects  of  furnace 
chamber  shape  on  natural  gas  in  an  industrial-scale  furnace  under 
MILD  operation.  The  furnace  chamber  shape  is  varied  by  changing 
the  angle  a  between  the  furnace  roof  and  sidewall  from  80°  to  100°. 
The  prediction  with  the  adopted  models  for  the  reference  case 
shows  good  agreement  with  the  experimental  results,  and  the 
conclusions  are  as  follows: 

In  the  present  furnace-burner  system,  the  furnace  chamber 
shape  has  a  significant  impact  on  in-furnace  mixing  and  combus¬ 
tion  characteristics.  Increasing  angle  a  would  cause  a  stronger 
recirculation  flow  field  with  more  combustion  products  entrained 
into  the  reactants.  This  results  in  a  broader  reaction  zone  and 
further  leads  to  a  lower  combustion  temperature  peak  and  NOx 
emission. 
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Fig.  12.  Contours  of  chemical  flame  shape  (a):  x-y  plane  (z  =  0);  (b):  x  =  1000  mm;  (c):  x  =  2000  mm.  The  black  curves  denote  the  level  of  Ro  =  0.99. 
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Fig.  13.  Chemical  flame  occupation  degree  for  the  five  cases. 


Table  4 

Predicted  NOx  in  the  flue  gas  for  the  considered  cases  (unit:  ppm,  dry). 


Casel 

Case2 

Case3 

Case4 

Case5 

a 

o 

o 

00 

85° 

o 

o 

o 

95° 

100° 

Total-NO* 

145 

140 

131 

126 

123 

Net-NO 

35 

30 

21 

16 

13 

Thermal-NO 

142 

138 

129 

125 

121 

Prompt-NO 

6 

6 

6 

7 

10 

Reburning-NO 

3 

4 

4 

6 

8 

The  Total-NO  is  calculated  by  activating  both  the  thermal-NO,  prompt-NO  and 
reburning  submodels.  While  Thermal-NO  and  Prompt-NO  are  calculated  separately. 


According  to  the  stability  limits  of  diffusion  combustion  pro¬ 
vided  by  Wiinning  et  al.  [4  ,  the  MILD  combustion  stability  in  the 
investigated  domain  is  found  to  be  enhanced  when  increasing 
angle  a ,  and  this  indicates  that  a  relatively  large  angle  a  is  helpful 
for  establishing  MILD  combustion  under  the  present  conditions.  In 
contrast,  the  combustion  regime  becomes  unsteady  when  angle  a  is 
less  than  85°,  and  this  further  implies  that  the  angle  a  cannot  be  too 
small  to  establish  MILD  combustion  under  the  present  conditions, 
or  the  initial  injection  conditions  must  be  intensified  enough. 
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